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Abstract 

The  micromechanics  involved  in  increased  crack  growth  resistance,  Kr,  of  a  HB2 
particulate-reinforced/SiC-matrix  composite  over  that  of  a  monolithic  SiC  was  analyzed 
experimentally  and  theoretically.  Fractography  studies  showed  that  the  advancing  crack 
was  attracted  to  the  particulates  due  to  the  tensile  region  surrounding  the  particulate. 
Countering  this  effect  was  the  compressive  thermal  residual  stress  in  the  matrix,  which 
resulted  in  the  toughening  of  the  composite.  This  thermal  residual  stress  field  in  a 
particulate  reinforced  ceramic  matrix  composite  is  induced  by  the  mismatch  in  the 
coefficients  of  thermal  expansions  of  the  matrix  and  the  particulate  when  the  composite 
was  cooled  from  the  processing  to  room  temperatures.  The  increase  'in,  Kr,  of  the  TiB2- 
particulate/SiC-matrix  over  the  monolithic  SiC  matrix,  measured  by  using  an  hybrid 
experimental-numerical  analysis,  was  77  percent,  and  compared  well  with  the 
analytically  predicted  increase  of  52  percent.  The  increase  in  Kr  predicted  by  the  crack 
deflection  model  was  14  percent.  Dependance  of  Kr  on  the  volume  fraction  of  the 
particulates  and  of  voids  is  also  discussed 
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I.  Introduction 

The  primary  toughening  mechanisms  of  particulate  reinforced  ceramics  have  been 
attributed  to:  (i)  interaction  between  the  crack  front  and  particulates  (crack  front  bowing 
model)1-2:  (ii)  crack  deflection  by  the  particulates  in  front  of  a  propagating  crack  (crack 
deflection  model);3  and  (iii)  crack  bridging  by  ductile  particulates  (particulate  bridging 
model).4  Other  secondary  mechanisms  which  contribute  to  the  toughening  of  ceramic 
composites  are:  (iv)  residual  stress  (strain)  field  due  to  mismatch  between  the 
coefficients  of  thermal  expansion  (CTEs)  of  the  ceramic  matrix  and  particulates;5-6  and 
(v)  the  grain  size.7 


Increase  in  the  fracture  toughnesses,  Kc,  of  particulate-reinforced/ceramic-matrix 
composites  with  modest  particulate  volume  fractions,  fp  =  0.1  -0.3,  over  the  fracture 
toughness  of  unreinforced  ceramic  matrix,  Km,  vary  from  10  percent7  and  40  to  70 
percent5-6-8-9.  This  increase  in  the  fracture  toughness  was  attributed  to  crack  front 
deflection.  Several  investigators  have  observed  such  crack  deflection  due  to  particulates 
in  particulate  reinforced  ceramic  matrix  composites.3-5-8  Such  crack  deflection  was  also 
observed  in  monolithic  ceramics  and  was  considered  a  major  contributor  to  the  toughening 
of  ceramics,10  although  the  crack  deflection  in  monolithic  ceramics  takes  place  along  the 
grain  boundary.  The  crack  deflection  observed  in  particulate-reinforced/ceramic-matrix 
composites,  however,  was  much  larger  in  magnitude,  and  was  thus  considered  a  dominant 
toughening  mechanism. 


The  crack  deflection  model  of  Faber  and  Evans3  is  based  on  a  geometrical  treatment  - 

of  a  crack  deflected  from  its  main  crack  plane  and  ignores  the  local  stress  field  at  and 
near  the  interface  between  the  matrix  and  the  filler  (particulate  in  this  paper).  This  local  >des 

stress  is  caused  by  the  difference  in  the  stiffnesses  and/or  the  coefficient  of  thermal  — - 
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expansion  (CTE)  between  the  matrix  and  particulate.  The  crack  deflection  model  predicts 
the  toughness  ratios  of  Kc/Km  *  1.12-1.15  for  uniformly  distributed  particles  and  for 
particle  volume  fractions  of  fp  -  0.1 -0.3.  This  K0/Km  ratio  predicted  by  the  crack 
deflection  model  is  smaller  than  1.4-1 .7,  typically  observed  in  most  particulate 
reinforced  ceramic  matrix  composites. 5,6,8, 9  Thus,  the  crack  deflection  model  only 
accounts  for  a  portion  of  the  increase  in  the  toughness.  An  additional  toughening 
mechanism  which  has  been  identified  in  fiber-reinforced/ceramic-matrix  composite  is 
fiber  bridging  in  the  region  behind  a  crack  tip.n-15  The  fiber  bridging  model,  however,  is 
only  applicable  to  short  and  continuous  fiber  composites  and  not  to  a  particulate 
reinforced  composite.  A  similar  idea  has  recently  been  used  by  Budiansky  et  al.  (1988)  to 
account  for  the  toughening  of  a  ceramic  by  the  bridging  of  ductile  metal  particulates, 
which  deform  elastically/plastically  until  failure.  This  model  is  valid  for  ductile  metal 
particulates,  but  is  not  valid  for  brittle  particulates  such  as  an  intermetallic  compound, 
for  example,  TiB2.  Moreover,  it  is  much  less  likely  that  TiB2  particulates  located  on  the 
crack  plane  will  break. 

An  additional  mechanism,  which  can  contribute  to  the  toughening  of  a  particulate- 
reinforced/ceramic-matrix  composite,  is  the  thermal  residual  stress  field  caused  by  the 
mismatch  of  CTEs  between  the  matrix  and  the  particulates.  The  role  of  thermal  residual 
stress  induced  by  CTE  mismatch  has  been  considered  as  a  major  cause  for  the  crack 
deflection  by  Wei  and  Becker,5  although  no  quantitative  estimate  on  the  increased 
toughness  was  made.  In  principle,  quantitative  analysis  of  toughening  due  to  crack 
deflection  by  thermal  residual  stress  can  be  made  within  the  framework  of  the  Fable* 
Evan's  crack  deflection  model,  but  such  an  analysis  would  become  intractable.  Another 
toughening  mechanism  that  has  not  been  studied  so  far  is  the  thermal  residual  stress  in 
the  matrix  domain.  If  such  residual  stress  in  the  matrix  is  to  contribute  to  toughening, 
then  the  thermal  residual  stress  in  the  matrix  must  be  compressive.  This  negative 
thermal  residual  stress  in  the  matrix  can  exist  when  the  CTE  of  the  particulate  exceeds 
that  of  the  matrix.  This  is  indeed  the  case  with  the  TiB2-particulate/SiC-matrix6-®.®  and 
TiC-particulate/SiC-matrix  composites.5 

In  order  to  quantify  the  influence  of  thermal  residual  stress,  a  combined 
experimental  and  theoretical  analysis  on  the  TiB2/SiC  composite  as  a  model  system  was 
undertaken.  Experimental  analysis,  followed  by  a  brief  description  of  the  analytical 
model  will  be  presented  in  Sections  II  and  III,  respectively.  Finally,  a  comparison 
between  the  experimental  and  theoretical  results  will  be  presented  in  Section  IV. 

ii.  Eapaiimflitt 

2.1  Materials 

The  ceramic  matrix  composite  considered  is  a  16  percent  volume  TiB2 
particulate/SiC  matrix  (TiB2/SiC)  composite.  The  unreinforced  SiC  and  TiB2/SiC 
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composite*  were  processed  by  pressureless  sintering  in  an  inert  atmosphere  at  a 
temperature  in  excess  of  2000°C.  The  room  temperature  material  properties  of  SiC 
matrix  and  HB2  particulate  are  given  in  Table  I.8-16-17 

The  temperature  dependency  of  the  CTE  of  SIC  and  TiB2  within  the  temperature 
range  of  200  to  2200K16-17  is  given  in  Fig.  1,  and  is  used  in  Section  III  to  calculate  the 
thermal  residual  stresses  induced  by  the  CTE  mismatch  between  the  matrix  and  the 
particulate. 

Photographs  of  typical  backscattered  electron  image  (BEI)  of  as-processed 
monolithic  SiC  and  TiB2/SiC  composite  are  shown  in  Fig.  2(a)  and  (b),  respectively.  In  Fig. 
2(b),  the  isolated  light  gray  areas  and  the  continuous  domain  of  gray  are  TiB2  particulates 
and  SiC  matrix,  respectively,  where  the  dark  isolated  spots  indicate  voids.  Voids  are  also 
seen  in  the  unreinforced  SiC  as  dark  spots  in  Fig.  2(a).  Some  of  the  voids  in  Fig.  2(b) 
appear  to  be  TiB2  particulates,  which  fell  out  during  polishing,  due  to  weak  bonding  of  the 
matrix-particulate.6  A  number  of  BEI  photos  similar  to  that  of  Fig.  2(b)  were  used  to 
measure  the  average  interparticulate  spacing,  X  (between  the  centers  of  two 
particulates),  and  the  volume  fraction  of  voids,  fv  (based  on  ASTM  Standard  E562-7618). 
The  results  are  X  *  11.1pm  and  fv  -  0.016,  which  will  be  used  in  the  analytical  model. 
The  value  of  fv  so  measured  is  consistent  with  the  density  of  the  composite,  i.e.,  (1— fv)  = 
0.98-0.99  reported  by  McMurty  et  al.6 

2.2  Measurement  of  Crack  Growth  Resistance,  Kr 

The  method  of  obtaining  crack  growth  resistance,  Kr,  is  based  on  a  hybrid 
experimental-numerical  analysis20  where  moire  interferometry  was  used.  The  purpose  of 
using  moire  interferometry  was  to  measure  the  crack  opening  displacement  (COO)  along 
the  crack  as  well  as  at  the  load  point.  A  schematic  of  the  moire  interferometry  system 
used  in  this  study  is  shown  in  Fig.  3.  A  moi  specimen  with  an  argon  laser  of  514.5  nm 
wavelength,  a  grating,  F,  of  600  lines/mm  ar.'„  ?  virtual  grating  of  2F  -  1200  lines/mm 
were  used.  The  displacement,  u,  can  be  measu.ud  by  counting  the  number  of  the  moire 
fringe  N19 


The  specimen  used  is  a  chevron-notched,  a  wedge-loaded  double-cantilever  beam 
(WL-DCB)  type,  which  is  shown  in  Fig.  4  where  the  dimensions  unit  are  in  mm.  A  side 
groove  with  a  depth  of  one  half  of  the  specimen  thickness  was  used  to  propagate  a 
straight  crack  along  the  crack  plane  (x-axis).  A  chevron  notch  was  used  to  generate  a  true 
crack.  The  effect  of  the  side  groove  on  the  stress  intensity  factor  is  minimal.21  The 
moire  interferometry  grating,  which  consists  of  a  thin  aluminum  film  on  an  epoxy 
adhesive,  is  attached  to  the  front  side  (not  the  groove  side).  The  displacement  along  the 


*  Supplied  by  British  Petroleum  Company,  Niagara  Falls,  N.Y. 
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y-axis,  Uy,  on  the  front  side  is  measured  by  the  moire  during  loading.  As  shown  in  Fig.  4, 
the  length  of  a  two-dimensional  crack,  a,  is  measured  from  the  intersection  of  the 
chevron  notch  with  the  front  side  (x  -  0,  y  -  0).  A  typical  moire  fringe  pattern  of  a 
TiB2/SiC  composite  specimen  is  shown  in  Fig.  S.  The  displacement,  uy,  along  the  x-axis, 
i.e.,  the  one-half  crack  opening  displacement  (COD)  was  measured  by  taking  the  average  of 
the  uy  values  on  the  left  and  right  side  fringes  since  the  fringe  pattern  is  often 
unsymmetrical  as  shown  in  Fig.  5 

For  the  hybrid  analysis,  the  uy  at  the  load  point,  i.e.,  the  contact  point  with  the 
WL-DCB  specimen,  together  with  the  measured  vertical  load,  px,  are  input  to  a  two- 
dimensional  finite  element  code  for  computing  the  stress  intensity  factor,  K|.  This 
procedure  eliminates  the  uncertainty  in  numerical  analysis  associated  with  the  friction 
at  the  loading  pin-specimen  contact.  The  boundary  conditions  and  the  pattern  of  the  finite 
element  analysis  are  shown  in  Fig.  6  where  only  one  half  of  the  WL-DCB  specimen  was 
used  for  the  analysis.  The  results  of  the  COD  measured  (symbols)  and  predicted  by  the 
finite  element  analysis  (solid  curves)  are  plotted  as  a  function  of  x  for  various  stages  of 
loading  in  Fig.  7.  This  finite  element  analysis  was  validated  by  the  excellent  agreement 
between  the  computed  and  measured  (by  moire  interferometry)  CODs  along  the  crack  plane 
in  Fig.  7.  Finite  element  analysis  was  then  used  to  compute  the  stress  intensity  factor 
using  the  standard  strain  energy  release  rate  method.  This  procedure  was  repeated  for 
each  incremental  increase  in  crack  length,  Aa,  which  then  yielded  the  crack  growth 
resistance  Kr. 

Crack  growth  resistance  curves,  KR-Aa,  were  thus  obtained  for  the  monolithic  SiC 
and  TiB2/SiC  composites  under  monotonically  increasing  load  as  shown  in  Fig.  8(a),  and 
cyclically  increasing  load  as  shown  in  Fig.  8(b).  The  resultant  KR-Aa  relations  are  shown 
in  Fig.  9  where  the  circles  and  the  squares  denote  the  monolithic  SiC  and  TiB2/SiC 
composites,  respectively.  The  open  and  filled  symbols  represent  the  results  of  monotonic 
and  cyclic  loadings,  respectively.  Fig.  9  shows  that  Kr  of  the  composite  increased  by  60 
percent  (cyclic  loading)  ~  80  percent  (monotonic  loading)  over  that  of  the  monolithic  SiC. 
Such  increases  in  Kr  cannot  be  predicted  by  the  existing  models  alone,  since  the  toughness 
increase  predicted  by  the  crack  deflection  model3  is  only  14  percent  above  the  monolithic 
matrix  and  is  much  lower  than  that  obtained  experimentally. 

iil  Analytical  .Model 

The  analytical  modeling  consists  of  two  steps.  The  first  step  is  to  calculate  the 
average  thermal  residual  stresses  in  the  matrix  and  reinforcement  by  using  the  Eshelby’s 
model22  which  is  modified  for  a  finite  volume  fraction  of  particulates.23-26  The  second 
step  is  to  calculate  the  change  in  the  mode  I  stress  intensity  factor,  AK,  due  to  the 
compressive  residual  stress  (obtained  in  the  first  step)  in  the  matrix.  In  this  model,  the 
crack  is  assumed  to  propagate  toward  the  nearest  tensile  stress  region,  which  coincides 
with  the  domain  surrounding  the  particulate  generated  by  the  thermal  residual  stress. 
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3.1  Thermal  Residual  Stress 


A  particulate  composite  is  assumed  to  consist  of  three  phases:  spherical 
particulates  of  volume  fraction  fp,  spherical  voids  of  volume  fraction  fv  and  the  matrix 
with  its  volume  fraction  of  1-fp-fv.  In  order  to  calculate  the  thermal  residual  stress  in  a 
composite  at  room  temperature,  Tr,  the  composite  is  assumed  to  be  subjected  to  a 
temperature  change,  AT  -  Tr-Tp  where  Tp  denotes  the  processing  temperature  such  as 
sintering  temperature.  The  thermal  residual  stresses  in  a  composite  are  assumed  to  be 
induced  by  elastic  deformations  of  the  matrix  and  particulates  under  a  uniform  AT.  The 
voids  are  included  since  the  SEM  observation  of  Fig.  2  shows  numerous  voids  distributed  in 
a  composite,  despite  the  small  volume  fraction  ~1.6  percent. 

The  modified  Eshelby's  model  is  suited  for  calculating  the  average  internal  stress 
field  in  a  composite.26-27  The  particulate  composite  containing  voids  subjected  to 
uniform  temperature  change,  AT,  is  schematically  shown  in  Fig.  10(a)  where  the  domains 
of  the  particulates,  voids,  and  the  entire  composite  are  denoted  by  £2j,  £22  and  D 
respectively,  thus  the  matrix  domain  becomes  D-£2i-£22.  The  elasticity  problem  of  Fig. 
10(a)  can  be  reduced  to  that  of  Fig.  10(b)  where  the  CTE's  and  aT  are  converted  to  the  CTE 
misfit  strains,  a1  (in  £2i)  and  g2  (in  £22).  which  are  given  by 


Tr 

g1  -  ^[(ap  -  am)8dT  (2) 

Tr 

g2  -  f(av  -  am)8  dT  (3) 

Tp 

and  where  ap,  av  and  am  are  the  CTEs  of  the  particulate,  void  and  matrix,  respectively,  and 
8  is  the  isotropic  tensor  (Kronecker's  delta).  The  subscripted  tilda  denotes  a  tensorial 

quantity.  av  is  not  involved  in  the  final  formulation  since  the  stiffness  of  £22  vanishes. 

The  modified  Eshelby's  model  Fig.  10(b),  is  equivalent  to  Fig.  10(c)  where  two 
inhomogeneities  (£2i,£22),  with  the  misfit  strains  defined  by  eqs.  (2)  and  (3)  are  replaced 
by  the  inclusions  with  eigenstrains  of  01  in  £2i  and  e2  in  £22.  The  thermal  stresses  in  a 
composite  can  then  be  solved  in  terms  of  the  eigenstrains.  Detailed  formulation  for  the 
thermal  stresses  in  a  particulate  composite  containing  voids  is  given  in  the  Appendix. 
Only  the  average  thermal  stresses,  which  will  be  used  later  in  this  section  for 
calculations  of  the  stress  intensity  factor  change,  AK,  will  be  considered.  The  average 
stress  field  in  the  particulate  <o>p  and  in  the  matrix  <o>m  are  isotropic  since  the 
particulate  reinforcement  and  the  CTE  misfit  strains  are  dilatational.  From  the  Appendix, 
these  are 


<o>p  -2(1-fp-fv)par 
Em  "  A 
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(5) 


2fpftai 
Em  "  A 


where 


(6) 


fp  and  fv  are  the  volume  fraction  of  the  particulates  and  voids,  respectively.  E|  and  vi  are 
the  Young's  modulus  and  Poisson's  ratio  of  Vie  i-th  phase,  where  i  «  m  and  p  represent  the 
matrix  and  particulate,  respectively.  When  the  composite  does  not  contain  voids,  the 
above  equations  can  be  simplified  as 


<q>m  2fppai* 
Em  "  A 


where 

A  -  (1-fp)(P+2)(1+Vm)  +  3pfp(1-vm)  (9) 

It  is  clear  from  the  above  equations  that  the  thermal  residual  stresses  are 
influenced  by  a  number  of  factors:  the  particulate-matrix  stiffness  ratio;  the  Poisson's 
ratios  of  the  matrix  and  particulate;  the  volume  fractions  of  the  particles  and  voids;  and 
the  CTE  misfit  strain,  (ap-am)AT.  The  formula  used  by  Wei  and  Becker  does  not  contain 
the  volume  fraction  of  the  particulates,  fp,  thus  it  is  a  lower  order  approximation  to  the 
stresses  in  a  composite  with  a  finite  volume  fraction  of  particulates. 

For  a  particulate  reinforced  ceramic  matrix  composite  with  ap  >  om,  such  as  a 
TiB2/SiC  composite,  ai*  <  0  as  seen  from  eqs.  (5)  and  (6)  the  average  thermal  stresses  in 
the  particulates  and  matrix  are  in  tension  and  compression,  respectively.  The  thermal 
residual  stress  field  in  a  TiB2/SiC  composite  consists  of  two  regions:  the  tensile  stress 
region  in  the  particulates  and  their  vicinity  in  the  matrix,  and  the  compressive  stress 
region  in  the  bulk  of  the  matrix,  as  shown  schematically  in  Fig.  1 1 . 

3.2  Stress  Intensity  Factor  Reduction  Due  to  Thermal  Residual  Stress  in  Matrix 

Consider  a  semi-infinite  crack  surrounded  by  a  particulate-reinforced/ceramic- 
matrix  composite  with  a  thermal  residual  stress  distribution  as  shown  in  Fig.  11.  The 
semi-infinite  crack  is  assumed  to  advance  from  one  concentrated  tensile  region  to  its 
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nearest  counterpart  which  may  be  off  the  crack  plane.  Thus,  the  crack  may  be  deflected 
according  to  Wei  and  Becher,6  but  the  general  crack  path  over  a  number  of  tensile  regions 
will  be  that  of  a  self-similar  crack  extension.  Hence,  a  self-similar  crack  extension  over 
an  average  Interparticulate  spacing,  X  will  be  considered.  When  a  semi-infinite  crack 
advances  by  a  distance  X,  which  is  the  Interparticulate  distance,  the  crack  tip  ligament  is 
subjected  to  the  local  average  compressive  stress,  q,  over  the  length  of  X-d,  as  shown  in 
Fig.  12,  in  addition  to  the  superposed  tensile  field  for  self-similar  crack  propagation.  "d" 
denotes  the  average  diameter  of  TiB2  particles.  The  existence  of  this  local  stress,  q, 
decreases  the  stress  intensity  factor,  AK|.  According  to  Tada,28  this  decrease  is 

AK,  «  |qV2*(X-d)  (10) 


IV.  Results  and  Discussion 
4.1  Results 

Using  the  analytical  model  developed  in  the  preceding  section,  the  increase  in  the 
fracture  resistance  of  a  TiB2  particulate-reinforced/SiC-matrix  composite  will  be 
predicted.  The  analysis  involves  two  steps:  calculation  of  the  average  thermal  residual 
stresses;  and  AK|. 

The  formulae  for  calculating  the  average  thermal  stresses  in  the  particle,  <o>p,  and 
the  matrix,  <o>m.  are  given  by  eqs.  (4)  and  (5),  or  eqs.  (7)  and  (8)  for  small  values  of  fv, 

respectively.  The  parameters  A,  (J  and  a,  must  be  calculated  first  using  known  input  fp, 

fv.  vm,  vp,  Ep  and  Em.  Since  the  measured  volume  fraction  of  voids  is  very  small,  i.e.,  fv  =* 
0.016,  eqs.  (7)  and  (8)  will  be  used  in  this  analysis. 

The  value  of  the  CTE  mismatch  strain,  a,,  is  defined  by  eq.  (2).  The  CTE’s  of  the 

matrix  (SiC),  am,  and  the  TiB2  particulate,  ap,  are  given  as  a  function  of  temperature  T(K) 
in  Fig.  1.  By  taking  Tr  -  300K  and  Tp  ■  2273K  and  integrating  op-am  with  respect  to  T 
over  the  entire  temperature  range  from  Tp  to  Tr,  eq.  (2)  yields 

aj  - -5.15  x  10-3  (11) 


Parameters  3  and  A  are  calculated  using  the  data  in  Table  1  and  are 
P  -  3.50 
A  -  6.86 


(12) 
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where  fv  «  0  was  assumed.  From  eqs.  (7),  (8),  (11)  and  (12)  and  using  fp  -  0.16,  the 
average  thermal  residual  stresses  are  computed  as 


<o>m  “  -345  MPa 
<o>p  -  1810  MPa 


(13) 


Next  an  estimate  of  the  toughness  increase  due  to  the  thermal  residual  stress  in  the 
matrix  is  made.  In  eq.  (10),  we  set 


q  -  <o>m  *  -345  MPa 


X  -11.1  pm 


(14) 


d  -  3.06  pm 

where  X  is  the  average  interparticulate  spacing  which  was  obtained  from  ten  SEM  photos 
containing  70  TiB2  particulates.  Substituting  eq.  (14)  into  (10)  yields 

AK|  ■  -1.56  MPa  (15) 

The  compressive  residual  stress  in  the  matrix  thus  reduces  the  stress  intensity  factor  by 
1.56  MPa.  This  reduction  in  AK|  is  equivalent  to  the  increase  in  the  crack  growth 
resistance  AKr  by  the  same  amount.  Hence 

AKr- 1.56  MPa  (16) 

4.2  Discussion 

The  measured  and  predicted  increases  in  fracture  resistance  is  now  compared.  The 
measured  increase  in  Kr,  AKr,  can  be  obtained  from  Fig.  9  for  the  case  of  monotonic 
loading, 


AKr  -  (Kr)  composite  -  (Kr)  monolith 

-  (5.3  -  3.0)  MPa  Vm  (17) 

«  2.3  MPa  VrrT 

It  is  clear  from  eqs.  (16)  and  (17)  that  the  major  portion  of  AKr  is  attributed  to 
toughening  due  to  the  thermal  residual  stress.  If  we  take  the  toughness  ratio  between  the 
composite, Kc,  and  the  monolithic  matrix,  Km,  the  measured  and  predicted  values  are 


Jk 

Km 


1  .77 
1.52 


measured 

predicted 


(18) 


The  measured  ratio  is  comparable  to  those  reported  by  others.6'8  The  crack  deflection 
model  of  Faber  and  Evans3  predicts  a  Kc/Km  ratio  of  1.14  for  uniformly  distributed  TiB2 
particulates.  Thus,  the  crack  deflection  model  cannot  account  fully  for  the  toughening  of 
a  ceramic  matrix  composite. 

The  predicted  value  of  K«/Km  ■  1.52  is  based  on  the  assumption  that  Tp  -  2273  K 
which  is  close  to  the  sintering  temperature  of  the  composite  and  also  that  both 
constituent  phases  of  SiC  and  TiB2  deform  elastically  during  the  cool  down  process  from 
2273K  to  300K.  SiC  a-phase  matrix  is  less  susceptible  to  creep  than  other  ceramics  at  a 
given  high  temperature  due  to  its  strong  covalent  bond  and  ionicity.  Its  activation  energy 
for  self  diffusion  is  large  and  results  in  a  slow  diffusion  rate  and  thus  slow  creep-rate.20 
This  slow  creep-rate  gives  rise  to  small  stress  relaxation.  It  is  thus  more  realistic  to 
assume  that  some  stress  relaxation  due  to  the  diffusion  takes  place  during  the  high 
temperature  regime.  If  T  -  1500°C  is  considered  as  a  threshold  temperature,  Tth.  above 
which  the  stress  relaxation  due  to  the  diffusion  can  occur,  then  Tp  must  be  set  equal  to 
1500°C  in  the  elastic  analysis  of  the  thermal  residual  stress.  For  Tp  -  1500°C,  the 
predicted  value  of  Kc/Km  -  1.43.  The  toughening  effect  based  on  Tp  -  1500°C  still 
constitutes  the  major  portion  of  the  observed  toughness  increase.  The  threshold 
temperature  Tth  is,  however,  difficult  to  determine  since  the  creep  (diffusion)  behavior  is 
a  continuous  function  of  temperature  and  time. 

Tth  of  metal  matrix  composites  can  be  estimated  by  in  situ  transmission  electron 
microscopy  (TEM)  study30  as  the  temperature,  at  which  the  dislocations  start  to  appear  in 
the  matrix  close  to  the  matrix-filler  interface  during  the  cool  down  process.  Such  an  in 
situ  TEM  study  of  a  ceramic  composite  is  difficult  to  conduct  since  Tth  is  too  high  for 
heating  stage  device  within  a  TEM.  In  addition,  the  generation  (punching)  of  the 
dislocations  from  the  matrix-filler  interface  due  to  the  CTE  mismatch31  is  unlikely  due 
to  high  friction  stress  of  the  ceramics  matrix. 

It  was  assumed  in  the  present  model  that  the  crack  propagation  path  over  a  number 
of  tensile  regions,  each  surrounding  a  particulate,  is  self-similar  despite  the  locally 
wavy  fracture  surface.  Figure  13(a)  shows  a  typical  BEI  photo  of  a  crack  propagating  in  a 
TiB2  particulate-reinforced/SiC-matrix  composite  toward  the  right  and  Fig.  13(b)  is  a 
portion  of  Fig.  13(a)  at  higher  magnification  .  It  is  noted  in  Fig.  13(b)  that  the  dark 
regions  on  the  crack  path  were  occupied  by  TiB2  particulates  which  fell  off  during 
polishing.  Both  figures  indicate  that  the  crack  actually  propagated  from  particulate  to 
particulate  located  on  the  crack  plane.  Figures  14(a)  and  (b)  are  BEI  photos  of  the  fracture 
surface  of  the  composite  and  monolithic  SiC,  respectively.  A  comparison  between  Figs. 
14(a)  and  (b)  reveals  that  the  degree  of  the  local  waviness  of  these  photos  is  the  same 
except  for  the  regions  next  to  particulates.  This  secondary  wavy  surface  which  is 
associated  with  the  particulates  are  probably  due  to  the  crack  deflection  proposed  by 
Faber  and  Evans3.  The  toughening  due  to  this  crack  deflection  was  evaluated  as  Kc/Km  - 
1.14.  By  superposing  the  toughening  due  to  crack  deflection  onto  that  due  to  compressive 
stress  In  the  matrix,  one  can  estimate  the  toughness  ratio  as 
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Km 


1  +  0.14  +  0.52  -  1.66 


(19) 


The  toughness  ratio  so  estimated  Is  close  to  the  experimental  value  of  Ko/Km  -  1.77. 

The  validity  of  the  present  model  can  be  partially  supported  by  the  experimental 
results  by  Jenkins  et  al.9  where  the  crack  growth  resistance,  Kr,  of  TiB2/SiC  composite 
measured  at  various  temperatures  from  room  to  1400°C  exhibited  a  decrease  with 
increasing  testing  temperature.  Toughening  in  the  present  model  is  proportional  to  the 
compressive  thermal  residual  stress  in  the  matrix,  which  in  turn  is  proportional  to  the 
temperature  change,  AT  ■  Tj-Tp,  where  Tt  is  the  testing  temperature  for  measuring  Kr. 
Hence,  the  higher  testing  temperatures  result  in  smaller  increase  in  Kr  as  observed  by 
Jenkins  et  al.° 

The  effect  of  volume  fraction  of  particulate,  fp,  can  be  examined  by  looking  at  eqs. 
(7),  (9)  and  (10).  The  thermal  residual  stress  in  the  matrix  <o>m  is  an  increasing  function 
of  fp,  as  shown  by  eqs.  (7)  and  (9),  and  will  influence  linearly  AK|  obtained  by  eq.  (10).  On 
the  other  hand,  the  average  interparticulate  spacing,  X,  is  proportional  to  l/Vfji31 

X-  1.085  d/VfJT  (20) 

while  AKi  is  proportional  to  Vx-d,  as  shown  by  eq.  (10).  Using  the  average  measured 
diameter  of  the  p articulates,  d  •  3.06  pm,  one  can  calculate  the  interparticulate  distance, 
X,  by  eq.  (20),  which  yields  X  -  8.14  pm.  This  value  of  X  is  smaller  than  the  average  of  the 
measured  interparticulate  distances  on  70  T1B2  particulates,  X  -  11.1  pm.  This 
discrepancy  between  the  measured  and  that  predicted  by  eq.  (20)  is  attributed  to  some 
error  involved  in  our  measurement  or  the  degree  of  approximation  inherent  in  eq.  (20). 
Non-uniform  distribution  of  TiB2  particulate  also  contributed  to  a  larger  measured  X  is 
larger  than  that  predicted  by  eq.  (20).  By  combining  the  above  effects,  the  toughness 
ratios,  Kc/Km,  are  computed  for  various  values  of  fp.  The  predicted  results  of  the  Kc/Km 
versus  fp  relation  are  shown  as  a  solid  curve  in  Fig.  15,  with  the  experimental  value  at  fp 
■  0.16  shown  as  an  open  circle.  Also  shown  in  Fig.  15  are  the  prediction  of  Kr  based  on  the 
crack  deflection  model  by  Faber  and  Evans  for  two  cases  of  particulates  distribution; 
uniform  spacing  (dash)  and  distributed  spacing  (dash-dot).  It  is  clearly  seen  from  Fig.  15 
that  the  present  model  can  predict  the  observed  toughness  increase  better  than  the  crack 
deflection  model  and  it  can  account  for  the  major  toughening  of  a  composite.  It  follows 
from  Fig.  15  that  the  crack  growth  rate  ratio  of  the  composite  to  the  monolithic  matrix 
increases  with  fp.  The  rate  of  increase  in  Kc/Km  is,  however,  not  linear  with  fp  but 
decreases  with  increase  in  fp.  This  fp  dependence  of  the  fracture  toughness  ratio  has  also 
been  observed  by  other  investigators.6’33  As  fp  increases  to  higher  values,  for  example, 
in  excess  of  0.5,  the  present  model  fails  to  predict  the  correct  trend  of  the  Kc/Km  versus 
fp  relation.  At  such  high  volume  fractions,  many  particulates  come  in  contact  and  form 
large  clusters.  The  microstructure  of  these  clusters  differs  from  that  of  the  present 
model,  as  shown  in  Fig.  10.  Thus,  the  Kc/Km  ratio  would  be  saturated,  or  even  start  to 
decrease  at  higher  volume  fractions.  Endo  et  al.33  have  indeed  observed  such  a  trend  with 
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several  non-oxide  particulates-reinforced/SiC-matrix  composite  systems.  Figure  16 
shows  the  Kic  versus  fp  relation  of  TiC  particulate/SiC  matrix  composite. 

The  present  model  predicts  high  tensile  stress  in  TiBz  particulate  at  room 
temperature  as  per  eq.  (13).  If  this  high  tensile  stress  <o>p  exceeds  the  bonding  strength 
of  the  matrix-particulate  interface,  ob.  then  the  interface  will  be  debonded.  It  seems 
appropriate  to  assume  that  the  interfacial  bonding  strength  varies  from  particulate  to 
particulate.  Thus,  some  particulates  with  ob  £  <o>p  would  become  debonded  voids  where 
the  debonded  interface  cannot  transfer  stress  except  for  compressive  stress.  It  is, 
however,  difficult  to  determine  the  percentage  of  such  particulates  with  weak  interface, 
but  one  can  obtain  a  rough  estimate  of  the  percentage  by  observing  the  percentage  of  TiB2 
particulates  that  would  fall  out  from  the  specimen  during  the  polishing.  If  the  debonded 
voids  so  nucleated  are  also  counted,  then  the  volume  fraction  of  voids,  fv,  increases  as  the 
fraction  of  TiB2  particulates  with  weak  interface  increases.  The  effect  of  fy  on  Kr  can  be 
studied  by  calculating  the  average  compressive  stress  in  the  matrix  <o>m  using  eqs.  (5) 
and  (6).  The  solid  curve  in  Fig.  17  shows  the  results  of  <o>m  as  a  function  of  fvby 
assuming  that  voids  are  nucleated  only  from  TiB2  particulates,  i.e.,  fv  +  fp  ■  0.16.  <o>m  is 
found  to  decrease  with  fv  and  becomes  zero  when  ail  the  TiB2  particulates  are  changed  to 
voids.  Thus  the  voids  nucleate  as  a  result  of  the  debonding  of  the  matrix-particulate 
interface  would  reduce  the  toughness  of  a  particulate  composite.  On  the  other  hand,  if  we 
assume  that  the  voids  are  induced  as  a  result  of  the  processing  (sintering  and  cool  down 
process),  and  the  matrix-particulate  interfaces  are  bonded,  then  fv  dependence  of  <o>m,  as 
shown  by  the  dashed  curve  in  Fig.  17,  has  the  opposite  trend  to  the  above  case.  Namely, 
<o>m  increases  gradually  with  the  increase  in  fv,  thus  contributing  to  the  toughness  of  a 
particulate  composite,  although  the  increase  in  the  toughness  is  modest.  The  values  of  fv, 
however,  are  usually  small  as  in  the  present  case  (fv  »  0.016).  Thus  the  effect  of  voids  on 
Kr  can  be  considered  minimal. 

The  effect  of  cyclic  loading  on  Kr  of  a  TiB2  particulate-reinforced/SiC-matrix  is 
obvious  from  Fig.  9.  Namely,  the  cyclic  loading  defined  by  Fig.  8  tends  to  reduce  the  value 
of  Kr  by  about  0.7  m-'/Pa.  This  reduction  in  Kr  can  be  explained  as  follows:  the  cyclic 
loading  modifies  the  local  stress  field  including  the  residual  thermal  stress,  resulting  in 
the  lower  compressive  stress  in  the  matrix.  Thus  a  crack  front  is  subjected  to  such  a 
modified  thermal  stress  field  during  the  loading  process  following  the  preceding 
unloading,  leading  to  the  reduction  in  the  toughness.  This  modified  thermal  stress  field 
would  not  occur  in  the  monolithic  SiC  due  to  the  lack  of  TiB2  particulates. 

iv.  Conclusion 

Increase  in  the  Kr  of  the  composite  over  the  unreinforced  matrix  is  attributed  to 
new  analytical  model  based  on  the  compressive  thermal  residual  stress  in  the  matrix  in  a 
particulate  reinforced  ceramic  matrix  composite.  The  Kr  predicted  by  this  model  agrees 
with  the  experimental  value  better  than  an  existing  model. 
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Table  1  Material  Properties  of  SiC  and  TIB2 


Unit 

SIC 

TIB2 

Crystal  structure 

hexagonal 

hexagonal 

Specific  gravity 

g/cm3 

3.21 

4.52 

Melting  point 

°C 

2540* 

2790 

Young's  modulus 

GPa 

410 

531 

Poisson's  ratio 

1 

0.19 

0.28 

Average  size 

pm 

grain  size  3 

2.06*** 

Volume  fraction 

0.84 

0.16 

CTE" 

1  o-«/oc 

4.02 

4.6 

*  Degassing  temperature 

**  These  CTE  values  taken  from  ref.  6  are  averages  over  room  temperature~700°C 
range,  and  are  not  accurate.  The  accurate  CTE  values  of  SiC  and  TiB2 
200~2200K16.17  are  given  in  Fig.  1. 

‘“This  value,  which  is  taken  from  ref.  6,  is  smaller  than  the  average  value  that  we 
have  measured.  This  is  discussed  later. 


14 


Appendix: 


The  Formulation  for  the  Thermal  Residual  Stresses 


Referring  to  Fig.  10(b),  the  internal  elastic  stress  field  in  domain  fii  (particulates) 
and  (i 2  (voids)  are  denoted  by  g1  and  g2,  respectively,  which  are: 


CP  •  (e  +01  - jj1) 

in  flii 

(A1 ) 

Cv  .  (e  +e2-a2) 

in  £22 

(A2) 

where  the  bold  face  letters  and  those  with  subscripted  tildas  denote  tensorial  quantities, 
CP  and  Cv  are  the  stiffness  tensor  (of  rank  4)  of  the  particulate  and  void,  respectively,  e1 

and  e2  are  the  strain  disturbed  by  Qi  and  £22,  respectively,  and  e  is  the  average  strain  in 
the  matrix,  and  related  to  the  average  stress  in  the  matrix  <o>m  as 

<g>m-Cm*e  (A3) 

The  dots  used  in  the  above  equations  denote  inner  product  between  tensors.  Following 
Eshelby,22  Fig.  10(b)  can  be  further  reduced  to  Fig.  10(c)  where  two  inhomogeneities, 
particulates  and  voids,  are  now  replaced  by  two  inclusions.  An  inclusion  is  a  domain 
surrounded  by  the  matrix  and  possesses  the  same  properties  as  the  matrix  and  also  an 
eigenstrain  e*  which  is  an  inelastic  strain.24  Thus  Equations  (A1)  and  (A2),  which  are 
formulated  based  on  Fig.  10(b),  are  now  reduced  to  the  following,  Fig.  10(c): 


gi  -  Cm  •  (e  +  e1  -  e*1)  in  Qi 

(A4) 

g2  «  cm  •  (e  +  e2  -  e*2)  in  (22 

(A5) 

Here  eqs.  (A4)  and  (A5)  are  valid  for  the  entire  composite  domain,  0,  if  the  eigenstrains, 
e*1  ana  e‘2,  are  defined  as  zero  in  the  domain  outside  fti  and  Q2,  respectively.  Noting 
that  the  integration  of  the  stress  field  g1  and  g2  over  0  vanishes: 

Q. 

< 

N 

O 

(A6) 

we 

obtain 

e--fp  (ei  -e*i)  -fv  (e2  -e*2) 

(A7) 
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where  fp  and  fv  are  the  volume  fraction  of  the  particulate  and  void,  respectively.  The 
strain  field,  e1  and  e2,  in  the  inclusions,  are  related  to  the  corresponding  eigenstrain,  a*1 
and  a*2,  respectively,  by  the  Eshelby's  tensor  S: 


a1  -  S  ’a*1  (A8) 

a2  -  S  •  a*2  (A9) 

Eshelby's  tensor  (of  rank  4)  for  an  isotropic  matrix  is  a  function  of  the  Poisson  ratio  of 
the  matrix,  vm.  and  the  geometry  of  fl.22-24 

By  equating  eqs.  (A1)  and  (A2)  to  eqs.  (A4)  and  (A5),  using  Cv  -  0  due  to  voids,  and 
eqs.  (A7)  through  (A9),  we  arrive  at: 

CP  •  l(l-fp)S  •  a*1  +  fpa*1  -  fv  (S  •  a*2  -  a*2)  -  q*1)  - 

Cm  •  I(i-fp)(S  -  a*1  -  a*1)  -  fv  (S  •  a*2  -  a*2)]  (Aio) 

-fpa*1  +  (1  -  fv)e*2  -  0  (Ail) 

where  the  isotropic  stiffness  tensors,  Cm  and  CP,  and  misfit  strain  q*1  are  given  in  index 
form  as: 


cj|ki  -  Xj8ij5ki  +  pj(8ik8ji  +  5||6kj),  J  «  m  or  p  (A12) 

*i  JR 

“IJ  -  J  («p  -  ®m)  8|j  dT  (A13) 

Tp 

where  8ij  is  the  Kronecker's  delta,  Xj  and  pj  are  Lam6  constants  of  the  J-th  phase 
material  (J  ■  m  for  the  matrix,  J  -  p  for  the  particulate),  am  and  ap  are  the  CTE  of  the 
matrix  and  particulate,  respectively,  and  Tr  and  Tp  are  the  room  and  processing 
temperatures,  respectively. 


The  problem,  which  is  illustrated  Fig.  10(c),  gives  rise  to  isotropic  state  since  the 
particulate  reinforcement  and  the  CTE  mismatch  strain  are  dilatatlonal.  Hence  only  one 
independent  component  for  each  unknown  eigenstrain  in  eqs.  (A10)  and  (All)  exists. 
Eshelby's  tensor  components  required  for  solving  eqs.  (A10)  and  (All),  thus  are: 


Sun 


S2222  -  S3333 


7  ~  5vm 
15(1  -  vm) 


Si  122  -  S2233  -  S331I-  Si  133  -  S3322  -  S2211 


5vm  -  1 
15(1  -  vm) 


(A14) 
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After  substituting  eq.  (A14)  into  eqs.  (A10)  and  (All)  and  eliminating  e*2,  we  can  solve 
*  *1  *1  *1 

for  e,  (-  e, !  -  %2  -  %3) 


(A15) 


where 


3Xp  +  2p.p  (1+Vm)  f'Ep"! 

2pm  (1-2vp)  (EmJ 


(A16) 


The  eigenstrain  in  n2,  a'2  (i.e. 


% 


*2 

®11 


*2 

%2 


can  be  obtained  from  eq.  (All)  once 


is  solved. 


The  average  stresses  in  the  particulate  <o>p  and  in  the  matrix  <o>m  are  related  to 
each  other  by  eq.  (A6). 

fp<o>p  +  (1-fp-fv)<o>m  -  0  (A1 7) 

where  the  stress  in  a  void  being  zero  was  used,  and  <o>j  with  J  -  m  (matrix)  or  p 
(particulate)  denote  an  isotropic  component,  not  tensor,  q 

The  average  stress  in  the  particulate  <o>p  can  be  calculated  from  eqs.  (A4)  and  (A7) 

•  tr 

through  (A9)  once  6|  and  efe  are  solved,  and  is  given  by 


<q>p 

Em 


-2(1-fp-fv)pa, 


<i-4{1'^Ki!)(1+,m) + 


(A18) 


The  average  stress  in  the  matrix  <o>m  titan  is  obtained  from  eqs.  (A17)  and  (A18) 

<o>m  _ _ 


Em  ...  Tf  fo_l._  ....  .  3Pfp(1-Vm) 

(1-fv)^|l-jf^|(P+2)(Uvm)  + 


(A19) 
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Eqs.  (A18)  and  (A19)  provide  the  average  thermal  residual  stresses  in  a  particle 
composite  with  voids  due  to  CTE  mismatch.  When  the  composite  does  not  contain  voids,  fv 
-  0,  eqs.  (A18)  and  (A19)  are  simplified  to 


<o>p 

Em 


Em 


-2(1-fp)Po,‘ 

[(1-fp)(P+2)(1+vm)  +  3pfp(1-vm)] 
2fppa1 

[(1-fp)(P+2)(1  +Vm)  +  3pfp(1-vm)] 


(A20) 


(A21 ) 
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Fig.  2 


Typical  backscattered  electron  image  of:  (a)  as-processed  SiC;  and  (b) 
TiB2/SiC  composite. 


DCB  Specimen  (Specimen 
—  Grating,  F  =  f/2) 


Fig.  3 


Moire  interferometry  setup 


loading  pin 


Wedge-loaded  double  cantilever  beam  (WL-DCB)  specimen  with  chevron  notch 
and  loading  pin.  All  dimensions  in  mm. 
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CRACK  GROWTH  RESISTANCE,  K„  (MPa-vftn) 
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TiB2/SiC  Composite 
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CRACK  EXTENSION,  A  a  (mm) 

Fig.  9  Crack  growth  resistance,  Kr  versus  crack  extension,  Aa,  relations  of  the 
monolithic  SiC  (squares)  and  TiB2/SiC  composite  (circles). 


Fig.  12  Semi-infinite  crack  advances,  X,  through  matrix  compressive  region  toward 
particulate  tensile  region. 


A 


Fig.  13  (a)  Backscattered  electron  image  (BEI)  of  typical  crack  propagating  path  in  a 

TiB2  particulate-reinforced/SiC-matrix  composite,  (b)  a  segment  of  (a)  at 
higher  magnification. 


Fig.  14 


BEI  of  fracture  surface:  (a)  TiB2/SiC  composite;  (b)  monolithic  SiC. 
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Volume  Fraction  of  Particulate  fp 

Fig.  15  The  crack  growth  resistance  ratio;  TiB2*particulate/SiC*matrix  composite 
to  the  monolithic  SiC,  as  a  function  of  the  volume  fraction  of  TiB2 
particulates,  fp.  Predictions  based  on  the  present  mode)  (solid  curve)  and 
Faber-Evans  crack  deflection  model  (dash  and  dash-dot  curves),  and  the 
present  experimental  result  at  fp  -  0.16  (open  circle). 
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Fig.  16  Kic  versus  fp  relation  of  TiC-particulate/SiC-matrix  composite.33 
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